function LPCVD_PP_ClosedLoop_r2_only(folder,r2_set,m2_set,repeat)
close all;

data_r2 = SinkToFile([folder,'/Output_r2_1.dat'],repeat);
data_r2.VariableName = '<r^2>';
data_r2.VariableUnit = 'layer^2';
figure(1);
plot(data_r2,'b-','linewidth',2);
if ~isempty(r2_set) 
    hold on;
    xlimit = xlim(gca);
    plot(xlimit,[r2_set,r2_set],'r:','LineWidth',2);
end
print('-djpeg',[folder,'/',folder,'_r2.jpg']);
print('-depsc2',[folder,'/',folder,'_r2.eps']);

data_m2 = SinkToFile([folder,'/Output_m2_1.dat'],repeat);
data_m2.VariableName = '<m^2>';
data_m2.VariableUnit = '';
figure(2);
plot(data_m2,'b-','linewidth',2);
if ~isempty(m2_set)
    hold on;
    xlimit = xlim(gca);
    plot(xlimit,[m2_set,m2_set],'r:','LineWidth',2);
end
print('-djpeg',[folder,'/',folder,'_m2.jpg']);
print('-dpesc2',[folder,'/',folder,'_m2.eps']);

data_W = SinkToFile([folder,'/Output_W_1.dat'],repeat);
data_W.VariableName = '<W>';
data_W.VariableUnit = 'layer/s';
figure(3);
plot(data_W);
print('-djpeg',[folder,'/',folder,'_W.jpg']);
print('-depsc2',[folder,'/',folder,'_W.eps']);

data_A = SinkToFile([folder,'/Output_A_1.dat'],repeat);
data_A.VariableName = '<A>';
data_A.VariableUnit = 'layer/s';
figure(4);
plot(data_A);
print('-djpeg',[folder,'/',folder,'_A.jpg']);
print('-depsc2',[folder,'/',folder,'_A.eps']);